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1  Introduction 


In  this  paper,  we  study  an  M/M/1  retrial  queue  with  an  unreliable  server  whose  orbit  and 
normal  queue  have  infinite  storage  capacity.  Customers  join  the  retrial  orbit  if  and  only  if  they  are 
interrupted  by  a  server  breakdown.  Retrial  customers  do  not  rejoin  the  normal  queue,  but  rather 
attempt  to  access  the  server  directly  at  random  intervals  independently  of  arrivals  or  other  retrial 
customers.  However,  these  interrupted  customers  can  regain  access  to  the  server  only  when  it  is 
operational  and  idle  and  repeat  service  until  they  have  been  successfully  processed. 

Our  primary  motivation  for  studying  this  model  stems  from  its  interesting  mathematical 
properties  and  its  distinction  from  other  commonly  analyzed  retrial  queues.  In  particular,  we 
are  able  to  analyze  the  interaction  between  the  orbit  and  the  normal  queue,  an  interaction  that 
does  not  exist  in  the  vast  majority  of  retrial  models  that  do  not  include  an  infinite-  or  nonzero- 
capacity  normal  queue.  Under  the  assumed  dynamics,  we  show  that  the  steady-state  orbit  size  and 
system  size  distributions  possess  a  stochastic  decomposition  property.  In  particular,  these  random 
variables  may  be  expressed  as  the  sum  of  two  independent  random  variables,  one  which  corresponds 
to  the  same  measure  for  an  instantaneous  feedback  model  (i.e. ,  one  with  an  infinite  retrial  rate) 
and  the  second  which  is  a  generalized  negative  binomial  random  variable.  Furthermore,  it  will 
be  shown  that  an  interesting  stability  result  emerges,  namely  that  the  normal  queue  may  remain 
stable  even  if  the  stability  condition  for  the  entire  system  (and,  in  particular,  for  the  orbit)  is 
violated.  Considered  here  are  two  types  of  breakdowns:  active  breakdowns  which  occur  during  a 
service  cycle,  and  idle  breakdowns  which  occur  while  the  server  is  non-failed  and  idle.  The  times 
between  customer  arrivals,  service  completions,  breakdowns,  retrials,  and  repairs  are  assumed  to 
be  exponentially  distributed  random  variables. 

Queueing  systems  with  breakdowns  have  been  studied  extensively  in  the  literature  as  have 
retrial  queues.  However,  the  literature  devoted  to  systems  with  both  retrial  queues  and  unreliable 
servers  is  comparatively  sparse.  The  seminal  papers  in  this  area  are  [1]  and  [11].  All  models 
considering  retrial  queues  subject  to  server  breakdowns  assume  an  M/G/l/1  loss  system  with 
the  exception  of  [6]  which  additionally  considered  a  finite-capacity  M/M /v/k  queue  for  which  the 
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author  derives  sufficient  ergodicity  conditions  but  does  not  provide  analytical  distributions  for  queue 
lengths  or  other  measures.  Although  [8]  considered  an  M/G/l  retrial  queue  with  infinite-capacity 
orbit  and  normal  queue,  the  authors  did  not  consider  a  server  that  is  subject  to  breakdowns.  In 
their  model,  arriving  customers  who  find  the  server  busy  may  join  either  the  orbit  or  the  normal 
queue. 

For  retrial  models  with  no  waiting  room  and  server  breakdowns,  customers  arriving  to  find 
the  server  unavailable  (busy  or  failed)  join  the  retrial  orbit.  Some  models  (cf.  [2], [3],  [18],  [13],  [14], 
[7],  [16],  [12]),  force  these  customers  into  the  orbit  while  others  ([11],  [6],  [5],  [4],  [9],  [17])  provide 
the  option  of  joining  the  orbit  or  departing  the  system.  With  the  exception  of  two  cases  ([3]  and 
[16]),  these  models  also  either  force,  or  provide  the  option  for,  in-service  customers  interrupted  by  a 
server  failure  to  join  the  orbit.  Our  model  is  distinct  in  that  arriving  customers  who  find  the  server 
busy  or  failed  join  the  normal  queue  whereas  interrupted  customers  always  join  the  orbit.  A  variety 
of  failure  types  are  considered  in  the  literature  including  starting  failures  ([18],  [13],  [12]),  vacations 
([7],  [14]),  active  breakdowns  ([6],  [16],  [17]),  and  like  our  model,  both  active  and  idle  breakdowns 
([11],  [2]  [3],  [5],  [4],  [9]).  Most  retrial  orbits  are  assumed  to  behave  as  infinite-server  queues  with 
identical  exponential  service  times;  however  some  models  (cf.  [7],  [12],  and  [17])  consider  orbits  as 
FCFS  queues. 

In  this  paper,  we  provide  the  steady-state  joint  distribution  of  the  orbit  size  and  normal 
queue  size  when  the  server  is  idle  (operational  and  non-occupied),  failed  (non-operational  and  being 
repaired),  or  busy  (operational  and  occupied).  From  these  results,  we  obtain  the  joint  probability 
generating  function  (p.g.f.)  of  the  orbit  size  and  normal  queue  size,  and  the  p.g.f.  of  the  overall 
system  size  (the  total  number  of  customers  in  orbit,  normal  queue  and  in  service),  independent  of 
the  server’s  status.  We  provide  a  necessary  and  sufficient  condition  for  stability  of  the  orbit  and 
system  and  a  necessary  condition  for  stability  of  the  normal  queue.  Moreover,  we  demonstrate 
the  stochastic  decomposability  of  the  orbit  and  system  size  distributions  and  provide  simple  and 
intuitive  expressions  for  the  limiting  distribution  of  the  server’s  status  and  standard  queueing 
performance  measures.  To  our  knowledge,  this  paper  is  the  first  to  present  such  results  for  a  retrial 
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queue  with  breakdowns  having  an  infinite-capacity  orbit  and  normal  queue. 

The  remainder  of  the  paper  is  organized  as  follows.  Section  2  provides  the  mathematical 
model  description  and  discusses  conditions  for  stability.  In  section  3,  we  provide  the  steady-state 
equations  and  the  limiting  distribution  of  server  status,  and  also  demonstrate  the  stochastic  de¬ 
composition  property  for  the  orbit  and  system  size  distributions.  Standard  queueing  performance 
measures  are  presented  in  section  4,  while  section  5  provides  illustrative  examples. 

2  Model  Description  and  Stability 


Customers  arrive  to  the  system  according  to  a  Poisson  process  with  rate  A  >  0,  and  service 
times  are  i.i.d.  exponential  random  variables  with  rate  ji  >  0.  Server  failures  occur  at  a  constant 
rate  £  >  0,  and  the  constant  rate  of  repair  is  a  >  0.  A  customer  interrupted  by  a  server  failure 
enters  the  retrial  orbit  and  spends  an  exponential  amount  of  time  there  with  rate  9  >  0,  after 
which  it  either  enters  service  (if  possible)  or  remains  in  the  orbit  for  an  additional  exponentially 
distributed  time  with  rate  6.  Denote  by  Qt  the  number  of  customers  in  the  normal  queue  at  time 
t,  and  let  Rt  denote  the  number  of  customers  in  the  orbit  at  time  t.  The  random  variable  Xt  is  the 
occupation  status  of  the  server  given  by 


Xt  = 


1,  if  the  server  is  busy  at  time  t 
0,  if  the  server  is  non-busy  at  time  t 


while  St  describes  the  operational  status  of  the  server  at  time  t  defined  by 

St  = 


1,  if  the  server  is  non-failed  at  time  t 
0,  if  the  server  is  failed  at  time  t 

The  continuous-time  stochastic  process,  {(Qt,  Xt,  Rt,  St)  :  t  >  0}  describes  the  state  of  the  system 

at  time  t.  Let  Nt  denote  the  total  number  of  customers  in  the  system  at  time  t  (i.e.,  in  orbit,  normal 

queue,  and  in  service).  The  process  {Nt  ■  t  >  0}  describes  the  evolution  of  the  system  size  over 

time.  Denote  by  {un  :  n  >  0}  the  Markov  chain  embedded  at  the  transition  epochs  of  {Nt  :  t  >  0}. 

We  assume  that  as  t  — >  oo,  Rt  =>  R,  Qt  =>  Q,  and  Nt  =>•  N,  where  (=$■)  denotes  convergence  in 

distribution. 

The  proportion  of  time  the  server  is  operational  is  a/ (a  +  £);  thus,  the  effective  service 
rate  is  an/ (a  +  £)  and  A  (a  +  Q/an  <  1  is  a  necessary  condition  for  system  stability.  As  we  shall 
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see,  this  condition  emerges  naturally  from  Corollary  1  of  section  3  and  Corollary  3  of  section  4. 
To  show  that  A(a  +  £)/a/r  <  1  is  also  sufficient,  we  need  to  prove  that  {vn  :  n  >  0}  is  ergodic 
when  A(a  +  ^)/ap  <  1.  It  is  easy  to  verify  that  {un  :  n  >  0}  is  irreducible  and  aperiodic;  thus  it 
remains  to  prove  that  it  is  positive  recurrent.  Pakes  [15]  proved  that  an  irreducible  and  aperiodic 
Markov  chain  {un  :  n  >  0}  is  positive  recurrent  if  |7fc|  <  oo  for  all  k  and  limsup^^  7&  <  0,  where 
7*.  =  E(vn+ 1  —  vn\vn  =  k).  In  our  model,  70  =  1  and  7*,  =  (A(a  +  £)  —  a/i)/(A(a  +  £)  +  ap),  for 
all  k  >  1.  Clearly,  if  A(a  +  £,)/ap  <  1,  then  |7fc|  <  00  for  all  k  and  lim sup^^  7^  <  0.  It  will  be 
shown  in  section  4  that  this  condition  is  more  restrictive  than  the  necessary  condition  for  stability 
of  the  normal  queue. 

Define  n k,i,j,i  as  the  limiting  probability  that  the  system  is  in  state  ( k,i,j,l ),  that  is 
™k,i,j,l  =  lim  P{Qt  =  k,Xt  =  i,  Rt  =  j,  St  =  l ) 

t — >00 

where  the  indices  k,  i ,  j,  and  l  correspond  to  the  normal  queue  size,  the  occupation  status  of  the 
server  (0  or  1),  the  orbit  size,  and  the  operational  status  of  the  server  (0  or  1),  respectively.  The 
transform  variables  27  and  z2  correspond  to  the  orbit  size  and  the  normal  queue  size,  respectively. 
Let 

OO 

$k,i,l(Zl)  =  ^  ^  *1  ^k,i,j,l 
3=0 

denote  the  generating  function  of  n with  respect  to  the  orbit  size  and  let 

00 

z2)  =  ^2  z\ 4>k, i,l(zi) 

k= 0 

denote  the  generating  function  of  <j>k,i,i(zi)  with  respect  to  the  normal  queue  size.  When  the  server 
is  idle,  failed,  or  busy,  we  respectively  denote  these  p.g.f.s  by  </>o,o,i(^i)>  V’o.oC-D,  z2),  and  i/ii^zi,  22 )• 
The  function  <^00,1(21)  is  the  first  derivative  of  ^>0,0, 1(^1)  with  respect  to  z\.  Define  as  the 

joint  probability  mass  function  (p.rn.f.)  of  R  and  Q  and  let  q(-)  denote  the  p.rn.f.  of  N.  Then 

OO  OO 

G(zi,z2)  =  k)ziz2  =  </>o, 0, 1(^1)  +  V’o, 0O1,  z2)  +  iIji,i(zi,z2) 

k= 0  j= 0 

is  the  joint  p.g.f.  of  the  orbit  and  normal  queue  size,  and 

OO 

H(z)  =  ^Z  qU)z3  =  <h,o,i(z)  +  ipoAz,  z)  +  ^1,1(2,  z) 

3=0 
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denotes  the  p.g.f.  of  the  system  size.  In  section  3,  we  derive  these  transforms  which  provide  the 
limiting  distribution  of  server  status  and  reveal  some  interesting  and  useful  stochastic  decomposition 
properties. 

3  Steady-State  Equations 

In  this  section,  we  derive  the  steady-state  joint  distribution  of  the  orbit  and  normal  queue 
size  when  the  server  is  idle,  failed,  or  busy,  respectively.  Subsequently,  we  demonstrate  the  stochas¬ 
tic  decomposability  of  these  distributions  independent  of  the  server’s  status. 


Theorem  1  The  generating  functions,  0o, 0,1(21),  ^0,0(21,  z2),  and  0i,i(zi,  z2),  are  given  by 


0o, 0,1(21)  =  ^  A(Q  +  ^P(Zl)c, 


(1) 


Vy  o(zi,z2)  =  p(zi)c+1 

x  ~£{A£(o  +  £)2i(zi  -  z2)  +  [an  -  A(a  +  £zi)]{(2i  -  z2)[£  +  A(1  -  z2)]  +  M1  -  z2)}} 

/u(a  +  £){[£z2  -  (n  ~  Az2)(l  -  22)][o  +  A(1  -  z2)\  -  o£zi} 

and 

ipi,i(zi,z2)  =  P(zi)c+l 

— A{(1  -  22)[o  +  £  +  A(1  -  z2)][an  -  A(o  +  £zi)]  +  (zi  -  22)£(o  +  £)[o  +  A(1  -  z2)]}  , 

li(a  +  £){[£22  -  (n  -  Xz2)(l  -  z2)][o  +  A(1  -  z2)\  -  o£zi} 

where  P(z\)  =  [an  —  A(a  +  £)]/[a/x  —  A(o  +  £zi)]  and  c  =  (o  +  £)/#. 

Proof.  Using  the  balance  equations  of  this  Markovian  system,  one  can  easily  obtain 


(a  +  A)00, 0,0(21)  =  £0o, 0,1(21)  +  £2100,1,1(21)  (4) 

(a  +  A)0fc, 0,0(21)  =  A0fe-!,  0,0(21)  +  £2i0fe,i,i(zi)  (5) 

(A  +  £)0o,  0,1(21)  +  6»zi0oiO1(2i)  =  a0o,o,  0(21)  +  At0o, 1,1(21)  (6) 

(A  +  n  +  £)0o,  1,1(21)  =  a0i,o,  0(21)  +  A00, 0,1(21)  +  //0i,  1,1(21)  +  00o,  0,1(21)  (7) 

(A  +  n  +  £)0fc,  1,1(21)  =  O0fc+1, 0,0(21)  +  A0fc_i, 1,1(21)  +  M0fc+I,i,i(2i)-  (8) 
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Multiplying  both  sides  of  Equation  (4)  by  z\  for  k  =  0  and  Equation  (5)  by  z\  for  k  >  1  and 
summing  over  all  k  >  0,  and  similarly  for  Equations  (7)  and  (8),  we  obtain,  respectively,  the 
following  two  equations: 

[a  +  A(1  -  22)]0O,o(2l,22)  =  £00, 0,l(2l)  +  £2i0i, 1(21,22)  (9) 

and 


[£22  -  (n  -  Az2)(1  -  22)]  01,1(21, 22)  +  [a0o, 0,0(21)  +  //0o,i, 1(21)]  = 

A2200, 0,1(21)  +a0o,o(2l,22)  +  022  0O;0, 1(21).  (10) 

Using  a  technique  employed  in  [6] ,  we  obtain  another  balance  equation  by  equating  the  flow  in  and 
out  of  the  set  En  =  {( k ,  i.  j.  1)  :  k  <  n  —  1},  n  >  1  which  leads  to 

M0n, l,l(2l)  +  a0n, 0,0(21)  =  A0n_i.  1,1(21)  +  A0n_ijOiO(2i),  U  >  1.  (11) 

Multiplying  both  sides  of  Equation  (11)  by  for  n  =  0  and  by  Z2  f°r  n  0  1,  summing  over  all  n, 
and  simplifying  we  obtain 


(M  -  Az2)0i.i(2i,  22)  +  (a  -  Az2)0o,o(2i,  22)  =  a0o, 0,0(21)  +  nfio,  1,1(21).  (12) 


We  first  obtain  an  expression  for  0o, 0,1(21)  by  setting  z^  =  1  in  equations  (6),  (9),  (10),  and  (12) 
and  solving 


,/  /  \  A£(a  +  £)  /  n 

The  general  solution  to  this  ordinary  differential  equation  is 


(13) 


0o, 0,1(21)  = - ,  (14) 

(a/j,  -  A(a  +  £21))  $ 

where  C  is  a  constant.  Now  substituting  (14)  into  equations  (9),  (12),  and  (13),  we  obtain  (1),  (2) 
and  (3)  up  to  the  multiplicative  constant  C.  The  final  results  are  obtained  after  the  normalization, 
0o, o,i(l)  +  0o, o(l,  1)  +  01, i(l,  1)  =  1,  which  leads  to 


C  =(ju(a  +  £))  1  [an  -  A(a  +  £)]c+1 
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wher  ec=  (a +  £)/#.  ■ 

Let  pi,  pf ,  and  ps  denote  the  limiting  probability  that  the  server  is  idle,  failed,  or  busy, 
respectively.  A  direct  consequence  of  Theorem  1  is  as  follows. 

Corollary  1  The  limiting  probabilities  pi,  pf,  and  ps  are  given  by 

Pi  =  lim  <£0,0,1  (zi)  =  -^7  - 
zi-*-i  a  +  £  p 

PF  =  lim  00,0(21,^2)  =  —7-7, 
zi~>l  a  +  4 

and 

pB  =  lim  0i  i(zi,  z2)  =  -• 
zi->l  p 

Z2~>1 

Corollary  1  reveals  that  A (a  +  £)/ap  <  1  is  necessary  to  ensure  system  stability  as  noted  in  section 
2.  Theorem  1  also  provides  the  means  by  which  to  obtain  the  joint  distribution  of  the  orbit  and 
normal  queue  size,  as  well  as  the  distribution  of  the  overall  system  size,  independent  of  server 
status. 

Corollary  2  The  probability  generating  functions  G(z\,z2)  and  H(z)  are  given  by 


G{z\,Z2)  = 


-[up  -  A(a  +  £zi)][a  +  £  +  A(1  -  z2)]{£(zi  -  z2)  +  (1  -  z2)[p  +  A(1  -  z2)}} 


+ 


p{a  +  £){[£22  -{p-  Az2)(l  -  z2)][a  +  A(1  -  z2)\  -  a£zi} 

_ — A£(a  +  g)(zi  -  z2)[a  +  £21  +  A(1  -  z2)] _ 

p(u  +  C){[^22  ~  (p~  Az2)(1  -  z2)][a  +  A(1  -  z2)\  -  a£zi} 


P{z  i)c+1  (15) 


and 


tt /  \  _  [ap  -  A(«  +  &Q][<*  +  £  +  A(1  -  z)\  .  ,c+1 

I1\z)  1  g\u..  \~\r~.  1  \/i  ~\i 


(16) 


(a  +  £){(p  -  Az)[a  +  A(1  -  z)]  -  A fzY 
Using  standard  methods,  (15)  and  (16)  can  be  used  to  obtain  the  mth  moment  (m  >  1)  of  R,  Q , 
and  N,  respectively,  as  well  as  their  limiting  distributions. 

The  stochastic  decomposition  property  has  been  observed  for  the  system  size  distribution  of 
many  M/G/l  models  including  those  with  vacations,  retrial  queues,  and  breakdowns  (cf.  [18],  [4], 
[9],  and  [12]).  The  property  implies  that  the  random  variable  of  interest  (e.g.,  orbit  or  system  size) 
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may  be  expressed  as  the  sum  of  two  independent  random  variables.  We  observe  that  Equations  (1), 
(2),  (3),  (15),  and  (16)  depend  on  the  retrial  rate  9  only  through  the  constant  c.  Allowing  9  — ►  oo 
yields  a  model  in  which  retrial  customers  instantaneously  attempt  to  re-access  the  server.  Conse¬ 
quently,  this  shows  that  the  orbit  (system)  size  is  the  sum  of  two  independent  random  variables: 
one  is  the  orbit  (system)  size  in  the  instantaneous  feedback  model  and  the  second  is  a  generalized 
negative  binomial  random  variable.  To  see  this,  note  that  each  of  the  aforementioned  expressions 
shares  a  common  multiplicative  factor  in  Theorem  1.  The  following  two  propositions  describe  the 
decomposability  of  the  orbit  and  system  size  distributions. 


Proposition  1  The  stochastic  decomposition  for  the  orbit  size  is  given  by 


where  r  =  A £/[a(n  —  A)]. 


(17) 


Equation  (17)  is  easily  verified  by  setting  Z2  =  1  in  (15)  and  simplifying.  The  left-most  term  on 
the  right-hand  side  of  (17)  is  the  generating  function  for  R  in  the  instantaneous  feedback  model 
(i.e.,  when  9  —>  oo),  and  the  right-most  term,  P{z\)c,  is  the  generating  function  of  a  generalized 
negative  binomial  distribution  ([10])  with  parameters  r  and  c.  Similar  behavior  may  be  observed 
for  the  system  size  distribution  in  the  steady  state  as  noted  in  Proposition  2 

Proposition  2  The  stochastic  decomposition  of  the  system  size  is  given  by 

H{z)  =  m  (rA^y '  <18) 

where 

=  [ap  -  X(a  +  gz)]  [a  +  £  +  A(1  -  z)\ 

(a  +  £)  {(/i  -  Xz)[a  +  A(1  -  z)\  -  A£z}’ 

The  quantity  B(z)  in  (18)  corresponds  to  the  generating  function  for  N  (the  overall  system  size) 
in  the  instantaneous  feedback  model.  In  section  4,  we  provide  standard  queueing  performance 
measures  for  this  system. 


4  Performance  Measures 


In  this  section,  we  use  (15)  and  (16)  to  obtain  the  standard  queueing  performance  measures 
noting  that  the  mean  system  size  and  sojourn  time  may  be  decomposed  into  three  components 
corresponding  to  the  server,  orbit,  and  normal  queue  measures. 


Corollary  3  The  limiting  mean  orbit  size,  normal  queue  size,  and  system  size  are  respectively 


given  by 


ptr>\  _  «A^[/r(/x  +  £  —  A)  +  A(a  +  £)] _ A £(a  +  £) 

n[ap  -  A  (a  +  £)]  [a(fj,  +  £)  —  A(a  +  £)]  6[ap  -  A  (a  +  £)]  ’ 


E(Q)  =  A 


+  £)  +  A(a  +  £)2 
p(a  +  £)  [ot(p,  +  £)  -  A(a  +  £)] 


F(AT\_  A[^£+(a  +  £)2]  A£(a  +  £)  ,  , 

1  ’  (a  +  £)[a/r- A(a  +  £)]  0[a/i  -  A(a  +  £)] ' 

Proof.  The  mean  orbit  size  is  obtained  by  evaluating  G'{z\ ,  1)  at  z\  =  1,  whereas  the  mean 
normal  queue  size  is  obtained  by  evaluating  G'(  1,  zf)  at  z 2  =  1.  In  a  similar  manner,  we  obtain  the 
mean  system  size  by  (16).  ■ 

As  £  — >  0  in  (19)  through  (21),  we  observe  that  E(R)  — >  0,  E(Q)  — >  A2///(/x  —  A),  and  E(N)  —* 
X/(n  —  A).  These  limiting  values  are  consistent  with  results  for  the  standard  M/M/1  queue.  More¬ 
over,  for  A,  /j,  a,  and  £  fixed,  the  mean  orbit  size  is  bounded  below  by 


mm  =  iim  E(m  =  aAy"+i„~/)+/<at/)1 

e^oo  —  \{a  +  £)][a(^  +  £)  —  A(a  +  £)] 


It  is  worth  noting  that,  if  retrial  customers  are  permitted  to  rejoin  the  normal  queue,  and  we  let 
6  — >  00,  the  model  converges  to  the  standard  M/M/1  queue  in  which  case  the  orbit  is  always  empty 
and  the  orbit  size  distribution  does  not  admit  a  stochastic  decomposition. 

By  inspection  of  Equations  (19)  and  (21)  we  observe  that  A(a  +  £)  <  apt  is  necessary  for  the 
stability  of  R  (and  N),  and  by  (20)  we  see  that  A(a  +  £)  <  a(/a  +  £)  is  necessary  for  the  stability 
of  Q.  That  is,  the  normal  queue  can  be  stable  even  if  the  orbit  size  stability  condition  is  violated. 
Owing  to  the  nature  of  the  orbit  dynamics,  retrial  customers  are  subordinate  to  normal  customers 
and  may  be  served  only  when  the  server  is  idle  and  operational.  Hence,  normal  queue  customers 


9 


experience  a  greater  effective  service  rate  than  do  retrial  customers,  and  thus,  it  is  possible  that 
the  orbit  may  continue  to  grow  while  the  normal  queue  remains  stable. 

We  further  note  that  the  steady-state  mean  system  size  E(N)  may  be  conveniently  decom¬ 
posed  into  its  constituent  elements,  namely  the  number  of  customers  in  orbit,  normal  queue,  and 
service.  The  mean  time  spent  in  orbit,  in  the  normal  queue,  and  in  the  system  are  obtained  by 
respectively  dividing  (19),  (20),  and  (21)  by  A. 

5  Numerical  Illustrations 


We  now  illustrate  the  behavior  of  the  mean  orbit  and  normal  queue  size  as  functions  of  the 
traffic  intensity,  failure  rate,  repair  rate,  and  retrial  rate.  In  Figures  1  and  2,  we  respectively  plot 
E(R)  and  E(Q ),  against  traffic  intensity  A//x  for  four  values  of  £  when  n  =  10,  a  =  2,  and  9  =  5. 
In  case  £  =  0,  the  model  is  equivalent  to  a  standard  M/M/1  queue  with  no  retrials  or  failures.  For 
£  =  0.5,  £  =  1.0,  and  £  =  1.5,  E(R)  rapidly  increases  to  the  saturation  point  of  a/(a  +  £).  Similarly 
in  Figure  2,  we  observe  that,  for  all  choices  of  £,  E(Q)  rapidly  increases  to  the  saturation  point  at 
which  a(^  +  £)  =  A(a  +  £)  (see  Equation  (20)).  It  is  obvious  that  the  orbit  becomes  unstable  more 
quickly  than  the  normal  queue  for  the  given  failure  and  repair  parameters. 


X 


X 

V- 


Figure  1:  Mean  orbit  size  (/x  =  10,  a  =  2,  d  =  5).  Figure  2:  Mean  normal  queue  size  (/x  =  10,  a  =  2,  9  =  5). 
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